Particle production at very low transverse momenta 
in Au+Au collisions at y/s NN = 200 GeV 
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We present results on charged particle production at very low transverse momenta in the 15% 
most central Au+Au collisions at ^/s NN — 200 GeV obtained with the PHOBOS detector at RHIC. 
The invariant yields were measured at mid-rapidity in the transverse momentum ranges from 30 to 
50 MeV/c for charged pions, 90 to 130 MeV/c for charged kaons and 140 to 210 MeV/c for protons 
and antiprotons. No significant enhancement in low transverse momentum particle production is 
observed as compared to extrapolations of identified particle spectra measured at an intermediate 
Pt range. The spectra tend to flatten at low pr, consistent with the expectations of transverse 
expansion of the system. 

PACS numbers: 25.75.-q 



Collisions of gold nuclei at the Relativistic Heavy Ion 
Collider (RHIC) provide the means to study strongly in- 
teracting matter under conditions of high temperature 
and energy density. The study of low-p^ particle pro- 
duction is particularly interesting as it is directly associ- 
ated with the long-distance scales, accessible in heavy ion 
collisions but out of reach in hadronic interactions. Any 
enhancement of the low-pr yields compared to extrap- 
olations from higher transverse momenta could indicate 
interesting effects, e.g. new long-wavelength physics phe- 
nomena jl|. Pion production may also be modified if a 
transient state with partially restored chiral symmetry is 
produced in the early stage of the collision [2, y, 3 . 

The spectra of identified hadrons measured at RHIC 
at t/s~^ = 200 GeV 0,0,0 tend to flatten in the low px 
region. This effect is more pronounced for heavier parti- 
cles and is commonly attributed to a collective transverse 
expansion of the system. Hydrodynamics-based models 
incorporating transverse expansion provide a satisfactory 
description of the spectra measured above 0.2 GeV/c 
0,0.01 Other mechanisms, such as initial state parton 
interactions [III fl2T ] or final state hadron reinteractions, 
can also describe the broadening of the spectra and an in- 
crease of (pt) with particle mass. Still, if transverse flow 
develops during the evolution of the system produced in 



heavy ion collisions, the effects of flow would likely be 
observed at very low pr- 

In this paper we present data on particle production at 
very low transverse momenta, pt < 0.2 GeV/c, a region 
thus far not explored at RHIC energies. The PHOBOS 
detector at RHIC has the unique capability to detect 
particles with very low px in the two-arm spectrometer 
13]. Each spectrometer arm consists of 16 layers of Si 
detectors. The six innermost layers are located in an ap- 
proximately field-free region near the beam pipe. The 
proximity of the sensitive detector layers to the interac- 
tion region, the small amount of material between the 
collision vertex and the Si layers, and the high segmen- 
tation of the Si detectors permit a precise measurement 
of particles with small pt ■ 

The data were taken during the 2001 RHIC run. The 
event triggering, determination of the collision centrality 
and the position of the collision vertex were as described 
in [lii. For this analysis we selected the 15% most 
central Au+Au collisions {(N part ) = 303, where N part 
is the number of participating nucleons). Events with 
the reconstructed vertex position along the beam axis 
between -7 cm and +14 cm, relative to the nominal in- 
teraction point were accepted. The centrality and vertex 
selections yielded 2 x 10 6 events. 
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In this analysis we searched for particles stopped in 
the 5 th spectrometer layer. The reconstruction of par- 
ticle trajectories and determination of the particle mass 
was based on a detailed analysis of the measured energy 
depositions in the first five spectrometer l ayer s. Selec- 
tion criteria based on the energy depositions [lj| for track 
candidates were derived from an analysis of single stop- 
ping particles obtained from GEANT simulations of the 
PHOBOS detector 16]. Initially, all possible 5-hit track 
candidates were formed from hits with large energy de- 
positions (dE/dx > 0.5 MeV normalized to 300/im of 
Si, roughly 6 times the deposition of a minimum ionizing 
particle). For the kinetic energies of particles stopping 
in the 5 th layer (about 8 MeV for pions, 19 MeV for 
kaons and 21 MeV for protons) the energy losses depend 
on the particle mass ^3 113 as shown in the inset in 
Fig. ^ Each track candidate satisfied the requirement 
\(dE/dx)i- < dE/dx >i | < 1.4(7;, i = 1,...,4, where 
the mean and the rms value <7;, are calculated from the 
[dE jdx)i distributions for single simulated stopping par- 
ticles. 

A mass parameter was defined as {M p )i = 
(dE/dx)i J2k=i Ek, for the first four layers, where Ek de- 
notes the energy deposited in the k th layer. For non- 
relativistic particles, the mass parameter {M p )i depends 
only on the particle mass, since (dE/dx)i ~ 1//3 2 and 
J2k=iEk ~ m/3 2 /2. Since the measured energy depo- 
sition in the fifth layer (£5) cannot be used, due to 
contributions from annihilation and decay products, we 
estimated £5 using the properties of the mass param- 
eter and the assumption that the particle stops in the 
5 th layer. Having estimated E 5 we calculate the aver- 
age mass parameter, M p , and the total deposited energy, 
Etot — Y^i=i^i- Fig- E shows a scatter plot of M p vs. 
E to t , for track candidates found in the data which satisfy 
the selection criteria on dE/dx in the first four Si layers. 

The momenta of stopped particles are restricted by 
imposing cuts on E to t (see Fig. QJ. For track candi- 
dates with emission angles smaller than 60°, not shown 
in Fig.^ about 10% higher upper limits on E to t were ap- 
plied as compared to those shown in Fig. ^ Additional 
cuts on the angular deviations of the candidate track 
from a straight-line trajectory were imposed to reject 
false tracks. The cuts are mass dependent to account for 
differences in the multiple scattering of low-momentum 
pions, kaons and protons. This analysis procedure was 
tested on samples of simulated low-momentum charged 
particles. Reconstructed single tracks simulated by 
GEANT were used to calibrate the final particle mo- 
menta. The transverse momentum resolution is esti- 
mated to be approximately 5%. These measurements are 
confined to the mid-rapidity region (—0.1 < y < 0.4) and 
cover a limited range of transverse momenta: 30 to 50 
MeV/c for charged pions, 90 to 130 MeV/c for charged 
kaons and 140 to 210 MeV/c for protons and antiprotons. 
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FIG. 1: (Mp,Etot) scatter plot for the candidate tracks af- 
ter the (dE/dx)i cuts in the first four spectrometer layers. 
Only track candidates with emission angles larger than 60° 
are shown. The upper limits on E to t are marked by vertical 
lines. An inset shows la bands on the specific energy loss in 
the first four Si layers for simulated particles. 



For further analysis the reconstructed particles were 
divided into narrow (Ay, Apx) bins. For each bin, cor- 
rections were applied to the raw yields to account for 
acceptance, reconstruction and particle identification in- 
efficiencies, and absorption in the beam pipe. These 
corrections were obtained by embedding simulated low- 
momentum particle tracks into real data events, and 
range between 0.2 and 0.7%. 

The acceptance and efficiency corrected particle yields 
are contaminated by various background contributions: 
feed-down from weak decays and contributions from sec- 
ondary, misidentified and ghost particles. For protons 
and antiprotons the dominant background contributions 
originate from secondary particles produced in the beam 
pipe or detector material and products from weak decays. 
To estimate the background contamination of the (p + p) 
yields, we analyzed the distributions of the measured dis- 
tance of closest approach {DC A) from the vertex posi- 
tion to the particle trajectory, determined by a straight 
line passing through the hits in the first and third spec- 
trometer layers. The experimental DC A distribution is 
broader and has a tail to larger DC A values compared to 
the same distribution for reconstructed primary proton 
tracks. Parameterized shapes of the DC A distributions 
for reconstructed primary, secondary and feed-down pro- 
tons, obtained from an analysis of HIJING events |19j . 
were used to describe the measured DC A spectrum. An 
example of the fit, assuming equal weights of the two 
background components, is shown in Fig. The fitted 
fraction of background particles, averaged over all possi- 
ble weights for the two background components amounts 
to 0.31 ± 0.04(stai.) ± 0A9(syst.). The large systematic 
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TABLE I: The estimated systematic errors (%) on the mea- 
sured invariant yields. 
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FIG. 2: DC A distribution for reconstructed proton tracks in 
the data sample (points with error bars). The result of the fit 
to the measured distribution (see text for more explanations) 
is denoted by the solid curve. The dashed (dashed-dotted) 
curve shows the contribution from primary (background) par- 
ticles. 



error assigned to this fraction reflects the dependence of 
the total background estimate on the relative importance 
of secondary and feed-down contributions. The fraction 
of background particles thus determined agrees with the 
estimate of 0.34 ±0.03, obtained directly from the recon- 
struction of HIJING events. 

To estimate the background corrections for kaons and 
pions we used HIJING simulations with yields tuned 
to account for the fact that HIJING overpredicts the 
yields of different particle species |20j. This adjustment 
of HIJING yields is important for evaluating the back- 
ground correction for charged kaons which is dominated 
by misidentified protons. The estimated background in 
the measured kaon sample is 0.10 ± 0.04, while the back- 
ground level in the unsealed HIJING events is higher 
(0.23 ± 0.03) due to the larger relative yield of protons 
to kaons. For protons and kaons the corrections are av- 
eraged over the full px range. 

The charged pion background estimates for scaled and 
unsealed HIJING events agree. In this case, feed-down 
from weak decays and ghost tracks are the main source 
of the background contamination. The large number of 
reconstructed pions in the simulated events permit the 
determination of the px dependence of the background 
corrections, which varies between 0.31 ± 0.04 and 0.22 ± 
0.02 for p T of 0.03 GeV/c to 0.05 GeV/c. 

The relative statistical uncertainties from the data as 
well as the statistical errors on the acceptance, efficiency 
and background corrections are 4 % for n ± , 10 % for K ± 
and 12 % for (p+p). Various sources of systematic errors 
were investigated, related to the backgrounds, the analy- 
sis method, and detector effects. Table H] summarizes the 
final estimates of the systematic uncertainties. 

The dominant contribution to the systematic errors is 
due to the uncertainty in the background estimates which 
are evaluated to be approximately 28% for (p + p) and 
about 10% for charged pions and kaons. 
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The uncertainties related to the analysis method were 
estimated by performing the reconstruction with varied 
selection criteria. There are also contributions estimated 
by the reconstruction of HIJING events, and of events 
which were generated with pt spectra approximately con- 
sistent with our measured data and with intermediat- px 
results from other experiments 0,0j0- In addition, a 
5% uncertainty in the pr scale and the sensitivity to the 
GEANT energy threshold contribute to the systematic 
uncertainty of the reconstruction method. 

The dominant contributions to the detector-related er- 
ror come from the differences in particle yields measured 
separately in the two spectrometer arms and under dif- 
ferent magnetic field polarities, and the calibration of 
energy depositions. Other detector effects such as mis- 
alignment of the spectrometer layers, and the accuracy 
of the vertex position, were found to be negligible. 

The individual contributions from different sources of 
systematic uncertainty are added in quadrature to derive 
the total systematic error for these measurements: 12% 
for n ± , 18% for if± and 32% for (p+p). 

The fully corrected invariant yields, 
(2ttpt)~ 1 d 2 N '/ 'dpTdy , of charged particles produced 
in the 15% most central Au+Au collisions at y/s NN — 
200 GeV are shown in Fig. [3| together with the yields of 
identified hadrons measured by PHENIX 6( for the same 
centrality selection. Different parameterizations have 
been used to extrapolate the measurements at higher 
Pt to pt ~ in order to determine the mean pt and 
integrated particle yields 0, IE • Measurements at very 
low pt provide some constraints on the choice of the 
functional form which fits the pt spectra. The functional 
form A[exp(niT /Tfi^ + e]^ 1 , withe = — 1(+1) for mesons 
(baryons), reasonably well describes the intermediate-pT 
PHENIX measurements (m T < 1 GeV/c 2 ). The fit 
parameter, T fit , equals 0.229 ± 0.005, 0.293 ± 0.01 and 
0.392 ±0.015 GeV/c 2 correspondingly for charged pions, 
kaons and (p + p). An extrapolation of the fits to low pt 
agrees with the measured 1ow-pt yields (see Fig. 0J. 

It was predicted that the tot spectra of different par- 
ticle species should follow a universal function of tot 
[11 , Il2l . One implication of the tot scaling described 
in 12j is that the spectral shapes should be independent 
of particle mass over a wide range of tot- In Fig. 01 the 
yields normalized at tot = 2 GeV/c 2 are shown as a 
function of tot- This choice of normalization was mo- 
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FIG. 3: Invariant yields as a function of pr- For 1ow-pt yields 
(closed symbols) the boxes show systematic uncertainties. For 
comparison the intermediate-pr measurements from PHENIX 
Q are depicted by open symbols. The fits to PHENIX mea- 
surements (solid curves) are extrapolated to low pr (dashed 
curves). See text for more details. 
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FIG. 4: Invariant yields, normalized at hit = 2GeV/c 2 , as 
a function of mr. Lower panel shows the tut dependence of 
the local inverse slopes. 



tivated by an analysis of the local inverse slopes, Ti oc , 
obtained from the exponential fits done with a sliding 
window enclosing 5 consecutive tot points. (The first 
inverse slope is from the fit to our average yield and the 
first four PHENIX points.) As illustrated in the lower 
panel of Fig. 0] these local inverse slopes converge at 
wit ~ 2 GeV/c 2 . The normalized spectra are con- 
sistent with the scaling hypothesis in the range of Tot 
from about 1.5 to 2.2 GeV/c 2 . However at low mr, close 
to the threshold, the spectra for heavier particles deviate 
from the charged pion spectrum, with larger deviations 
observed for the (p+p) yields than for the charged kaons. 



This observation is consistent with mass dependent flat- 
tening of the spectra at low tot, expected in the presence 
of transverse expansion of the system, and contradicts the 
suggested tot scaling [l2j |. 

In summary, the yields of identified particles have been 
measured at very low px (below 210 MeV/c) in central 
y / s^" = 200 GeV Au+Au collisions. We see no enhance- 
ment in the production of particles with very low trans- 
verse momenta, which might have indicated the pres- 
ence of unusual long-wavelength phenomena. The low- 
Pt yields agree with extrapolations from intermediate-pT 
measurements. Finally, the measurements in this region 
of phase space, unexplored until now, show clear flatten- 
ing of particle spectra that increases with particle mass. 
This observation supports an interpretation of the data in 
terms of a collective transverse expansion of the system. 
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